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Nuclear effective interactions are often modelled by simple analytical expressions such as the 
Skyrme zero-range force. This effective interaction depends on a limited number of parameters that 
are usually fitted using experimental data obtained from doubly magic nuclei. It was recently shown 
that many Skyrme functionals lead to the appearance of instabilities, in particular when symmetries 
are broken, for example unphysical polarization of odd-even or rotating nuclei. In this article, we 
show how the formalism of the linear response in infinite nuclear matter can be used to predict and 
avoid the regions of parameters that are responsible for these unphysical instabilities. 



I. INTRODUCTION 
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' Methods based on the mean-field approach represent the tool of choice for the study of all atomic nuclei but the 
QQ lightest. A fairly good approximation of the mass and size of doubly magic heavy nuclei can be obtained at the mean- 
Ps| [ field level, often called the Hartree-Fock (HF) approximation, while the properties of other nuclei can be obtained 
by using a hierarchy of refinements from the inclusion of pairing correlations to the description of the fluctuations of 
' f— I ' ' collective degrees of freedom using the Generator Coordinate Method (GCM) and projection techniques (for symmetry 

, restoration). 

_Jh ' The key object to build is the Energy Density Functional (EDF). Usually, its expression is derived from an effective 
^ [ interaction like the GognyfT], 0| or Skyrmejl, [j] ones that depend on a small number of parameters. The Skyrme 
interaction being of zero-range type, the obtained EDF depends on local densities only when the usual Slater approx- 
imation is adopted for the Coulomb contribution. The calculation of the average value of the Skyrme interaction in a 
vacuum that is a product of occupied quasiparticle states generates all the terms of the functional in a consistent way, 
i.e. the contributions of the normal and anomalous (pairing) densities as well as the contributions of the time-even 
and time-odd densities. 

On the one hand, the effective character of this approach justifies the omission of certain terms in the functional (as, 
^.f^ ' for example, the so called "J^" terms) or the use of different parameters for its normal and pairing parts. Moreover, 
\^ since most of the time, the procedures used to adjust the parameters of the Skyrme interaction use only informations 
from time-even systems, one can also neglect several time-odd terms in the functional. When such a choice is done, 
the backward link between the functional and an effective interaction is broken. On the other hand, it was recently 
reahzed that the use of a functional not exactly derived from an effective interaction leads to spurious contributions 
in the energy in calculations beyond the mean-field approximation. These spurious contributions can be removed 
using a procedure recently formulated. [5, ^] It was also shown that if the effective interaction does depend on non 
integer powers of the density, other spurious contributions pollutes the energy in beyond mean-field calculations and 
the aforementioned correction procedure can not be used.[7| 
I The use of a single effective interaction with only integer powers of the density (or with explicit three and four body 
interactions) to generate all the terms of the functional would be the ideal choice to avoid or to be able to handle 
correctly the pathologies in beyond mean-field calculations. Such interactions do not exist at this moment and their 
developments are only at their early stages. 

The present work is devoted to the study of finite size spin and isospin instabilities in nuclei that appear at the 
mean-field level. We use a set of effective interactions that are built following the same protocol as the one used 
for the famous SLy interactions. 10] These interactions are supplemented by an effective two-body zero range tensor 
term. Since only doubly magic nuclei and infinite nuclear matter in its normal phase are calculated, the choice for the 
interaction in the pairing channel is not in question here. The parts of the functional that depend on the time-odd 
densities are obtained from the effective interaction and are all kept in the calculations. Having a non-integer power 
of the isoscalar density Pq in the density dependent term is not a problem since we do not perform beyond-mean-field 
calculations here. 



* Based on talk presented at 18th Nuclear Physics Workshop "Maria and Pierre Curie", 2011, Kazimierz, Poland. 



2 



This article is organized as follows. In section |TT] we present and discuss the form of the Skyrme interaction used 
in this study and its parameters fitting procedure is recalled in section IIIII Examples of finite size instabilities are 
shown in section ITVl Finally, in section |Vl the linear response formalism is briefly exposed and its ability to predict 
and avoid the regions in the parameters space that lead to instabilities is demonstrated. 



II. THE SKYRME EFFECTIVE INTERACTION 



The interactions and functionals used in this article to illustrate the appearance of finite size instabilities are: ( i) a 
serie of functionals obtained by modifying only one coefficient on the functional derived from the SLy5 interaction [lOj 
and leaving the other coefficients unchanged; (ii) the serie of TIJ interactions built for the study of the effect of tensor 
couplings on the properties of even-even nuclei. prilll2j 

In the TIJ intaractions, the central and spin-orbit parts are of the same form as for the SLy interactions, the tensor 
part being 

VT(r) ^\t, {5{y) [3 ((Ti . k) {CT2 ■ k) - ((Ti . rT2) k2] + h.c.} 

+ \t,{[i (O-I • k') 5{V) {CT2 ■ k) - (*Tl • rT2) k' • J(r)k] + h.c.} , (1) 

where h.c. stands for "hermitian conjugate". The parameters and along with ii, xi, t2 and X2 determine 
the eight coupling constants C^, Cf , C^'^ and C^'^ (with t = and 1 for the isoscalar and isovector parts of the 
functional). The definitions of the operators k, k', cri and (T2 as well as the list and expressions of the coupling 
constants appearing in the functional can be found in the article of Lesinski et al. [ll[ 



III. PARAMETERS FITTING PROCEDURE 



The method used to fit the parameters of a Skyrme interaction or functional consists in constructing a merit 
function that measures the deviation between targeted and calculated quantities. These quantities are generally 
several properties of the infinite nuclear matter, masses and charge radii of a set of doubly magic nuclei and one 
splitting between two single particle states, this latter being added to optimize the spin-orbit strength. 

For the serie of functionals built from the SLy5 interaction, we have modified the coefficient C^'' given by 

^^f^ + + + . (2) 

With the parameters of SLy5 the value of this coefficent is 16.375 MeVfm^. The serie of functional used here is then 
obtained by varying this coefficient without modifying the other ones. Changing a single coupling constant without 
refitting the entire functional deteriorates the properties of the functional itself {i.e. the reproduction of some physical 
observables) , but this is not crucial in the present work since this modification is only used to explore the relation 
between the appearance of instabilities and the value of one term of the functional. 

The set of T/ J interactions are obtained by setting = 30( J + / - 4) MeV fm^ and = 30( J - /) McV fm^ for 
/ and J integers from 1 to 6 and with 

Ccf = -1*1 {^1 + 5) + I *2 (a;i + i) - i {te + 3t„) , (3) 

= -J^h + j^t2- \ {te- to) . (4) 

AH the parameters of the interactions being then adjusted in order to minimize a given merit function [TT| with an 
iterative procedure nearly identical to that used for the construction of the SLy parameterizations. [lOj 



IV. INSTABILITIES IN FINITE NUCLEI 



To adjust the parameters of an effective interaction, one starts from an initial set of parameters that can at least 
provide with a saturation point in symmetric infinite nuclear matter and this set of parameters is modified iteratively 
in order to minimize a merit function. During the iterations, some instabilities can be experienced. For example, 
one can observe domains in finite nuclei where the isovector density becomes large. This phenomenon can ultimately 
lead to a total separation between protons and neutrons and a disappearance of nuclei as bound states. The regions 
of the parameters space that lead to these unwanted situations should be avoided and it would be helpful to have a 
criterium that predicts where those regions are. 
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In a previous work, [l3j it was shown that the coefHcient of the Skyrme functional can be related with this kind 
of isospin instability. We illustrate this correspondance by the following schematic serie of calculations for ""^Ca. We 
start from the effective interaction SLy5 which is stable in the isospin channel and for which C^'' — 16.375 MeVfrn'"^. 
We keep all the other coefficients of the functional constant and vary C^'' from 16 MeVfm^ to 37.1 MeVfm^ where 
an isospin instability appears in this nucleus (the limit value can change from one nucleus to an other). 

We observe that the total energy of ^®Ca calculated at the HF approximation decreases with C^'' (from -415.9 MeV 
to -An .2 MeV) and that this gain of energy is due to the appearance of strong oscillations of the isovector density 
Pi — Pn — Pp as shown on Fig. [TJ For C^'' > 37.1 MeVfm^ the nucleus ''^Ca is disrupted in finite size domains of 
neutrons and protons. 

This kind of instabilities can be detected during the fit procedure but the situation is obviously more complicated 
for the instabilities in the spin channel since one has to break the time reversal symmetry to observe it. A possible 
constraint on the Landau par ameters (as it was the case for the SLy functionals) would only avoid global polarizations 
of infinite nuclear matter, [ij] but not the finite size instabilities. 

Unfortunately, these spin instabilities have been identified recently in mean-field calculations where the time reversal 
symmetry is broken (cranked calculation [l^ and calculations of odd-even nuclei [ll|) when the functional is explicitely 
derived from the effective interaction {i.e. all time-even and time-odd terms are kept). Calculations of systems where 
the time reversal symmetry is broken are two much time consuming to be easily incorporated in the fit procedure of 
the parameters, one needs an efficient tool to predict and avoid all possible finite size instabilities both in the spin 
and isospin channels. 



V. LINEAR RESPONSE 



A. Formalism 



The basic formalism of linear response theory has been already explained in several articles, in our work we follow 
the formalism used in the article of Garcia- Recio et a?.[l3] In this work, we restrict our study to the response to 
an external probe of an infinite system initially unpolarized both in spin and isospin spaces (symmetric nuclear 
matter, SNM) at zero temperature. The excitations of the system result from the residual particle-hole interaction 
between particles below and above the Fermi level. In the theory of Fermi liquids, this interaction is given by the 
second functional derivative of the Energy Functional E[p], with respect to the densities taken from the Hartree-Fock 
solution. Its matrix elements read 

P" (Sp(ki -I- q, ki )(5p(k2 , ka + q) 

where we denote by a = (S,M;I) the spin and isospin in the particle-hole (p-h) channels with S = (or 1) for the 
non spin-flip (or spin-flip) channel, / = (or 1) for the isoscalar (or isovector) channel, and M being the projection 
of the spin on the quantization axis. As usual, we have chosen as independent variables the initial (final) momentum 
ki (k2) of the hole and the external momentum transfer q. To calculate the response of the medium, it is convenient 
to introduce the Green's function, or p-h propagator, at the energy w: G^"' (q, w, ki), which in the Hartree-Fock 
approximation, does not depend on the spin- isospin channel (a) and reads [Tsj 

( V w g(fcf -fci)-g(fcf -|ki + q|) 
uj + e{ki) - e{\ki + q|) + zyyo; 

Considering now the residual interaction of equation ([5]), we can calculate the RPA correlated Green's function by 
solving the Bethe-Salpeter equation 

GnPAili^^'^i) = GHF('Z,w,ki) 



(7) 



For a detailed discussion on the solution of equation ([7]) we refer to the work of Davesne et al. [l9| Finally the response 
function Xrpa(<Zi'^) of the infinite medium is 

d^fci 



(8) 
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where g is the spin-isospin degeneracy factor ((7 = 4 for SNM). Within this formahsm, an unstable mode is marked 
by the occurrence of a pole in Xr^a at o; = 0, corresponding to zero excitation energy. The occurrence of a pole at 
finite q characterizes a system that is unstable within respect to the appearance of a spatial oscillation of a given type 
with a given wavelength A ^ Stt/^/. 



B. Link with finite size instabilities in nuclei 



We search for the density pcrit in the (q, po) plan that corresponds to a pole in equation ([5]) at zero excitation energy. 
These densities are represented on the right panel of Fig. [T] for three modified SLy5 functionals in the 5* = 0, / = 1 
channels. 




FIG. 1: Left panel: Isovector density pi for Cf" = 16 MeVfm^ and Cf" = 37.1 MeVfm'"^ in **Ca. Right panel; Critical 
densities for three versions of the modified SLy5 functionals (Cf " = 16, 37 and 47 MeVfm^) in the S = 0, / = 1 channel. The 
horizontal dotted line represents the saturation density of the system paat- 

One can see that there is a clear relation between the instabilities discussed in section ITVl and the critical densities 
represented on the right panel of Fig. [TJ in particular we observe that the poles of the response function in the scalar- 
isoscalar channel approach the saturation density psat = 0.16 fm~^ for increasing values of C^'' . It is important to 
remark that this particular instability can not be detected using the Landau parameters since they only probe the 
limit 5 — > 0. 

The instabilities in the spin channels can not be shown as easily as the isospin ones since one has to break the time 
reversal symmetry and therefore the spherical symmetry as well. In Fig. [2] we show the minimal value pmin of the 
critical densities pcrit for any q for all the T/J interactions in the channels with S — \. To have a stable interaction, it 
is necessary that the functional satisfies the minimal condition pmin > Psat in all channels. As a result we observe that 
none of the T/J interactions is stable, this is in agreement with recent cranked HFB calculations in finite nuclei[15]. 



VI. CONCLUSION 



We have given and illustrated an example of the finite size spin and isospin instabilities that can occur in nuclei. 
We have also reported that finite size instabilities are known to appear in calculations that break the time reversal 
symmetries with several Skyrme interactions. 

The terms of the functional that are responsible for the appearance of the instabilities can be indentified and omitted, 
but this rough method breaks the link between the functional and an underlying effective interaction. Calculations 
beyond the mean-field approximation are then more involved since one has to remove spurious contributions from the 
energy. If one wants to keep the link that relates the interaction to the functional and avoid those spurious contributions 
in the energy, one needs a method that can predict there appearance and this method must be numerically fast enough 
to be usable during the fitting procedure of the parameters. 

The method of the linear reponse that we present here seems to be an efficient tool to reveal these finite size 
instabilities, although the exact correspondance between the instabilities in infinite nuclear matter and nuclei needs 
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FIG. 2: Minimal value of the critical densities in the four channels with 5 = 1 for the TIJ interactions. The correspondance 
between the symbols and the values of I is given on the top left panel. The J index of the interactions corresponds to the 
coordinate on the horizontal axis. The horizontal dotted lines represent the saturation density of the system psat . 



to be studied in more details. First of all, finite size and shell effects being non trivial in nuclei, one has to test whether 
there is a robust correspondance between instabilities in infinite nuclear matter and nuclei. If this is confirmed, one 
has to estimate how close the critical densities can be from the saturation density before an instability shows up. 

The calculations of the linear response and critical densities being very fast, they will be incorporated in the fitting 
procedure of effective forces and functionals in the future. 
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